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This  report  is  based  on  five  lectures  given  by  Dr.  I.  N.  Sneddon  at 
North  Carolina  State  in  April,  1963.  The  research  reported  here  concerns 
certain  crack  problems  in  the  mathematical  theory  of  elasticity.  A  group 
of  these  problems  were  presented  to  us  by  Dr.  George  Irwin  of  NRL  at  the 
beginning  of  this  sponsored  research  work;  the  present  report  concerns 
the  solution  of  some  of  these.  Part  of  the  research  results  presented  in 
this  report  have  already  been  submitted  in  the  form  of  three  papers  for 
possible  publication. 

Copies  of  this  report  are  being  distributed  as  directed.  This  project 
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THE  EFFECT  OF  INTERNAL  CRACKS  CN  THE  DISTRIBUTION  OF  STRESS  IN 
THIN  ELASTIC  STRIPS  AND  CYLINDERS 


1 .  Introduction. 

The  present  report  is  based  on  a  series  of  five  lectures  given 
in  the  Mathematics  Department  of  North  Carolina  State  College  in 
April  1963*  in  which  there  was  presented  a  connected  account  of  some 
recent  researches  in  the  classical  (infinitesimal)  theory  of 
elasticity.  In  particular  the  work  done  by  the  author  and  three 
collaborators  R.  P.  Srivastav,  R,  J,  Tait  and  J.  L.  Welch  in  the 
University  of  Glasgow  during  the  preceding  year  is  discussed  in  some 
detail. 

In  §2  there  is  a  discussion  of  the  two-dimensional  problem  of 
determining  distribution  of  stress  in  a  very  long  strip  of  uniform 
width  2  c  which  has  a  Griffith  crack  of  length  2  a  situated  symmetrically 
in  its  interior.  The  stress  field  is  supposed  to  be  set  up  by  the 
application  of  a  known  pressure  to  the  inner  surfaces  of  the  crack. 

The  equations  of  plane  strain  are  used  throughout  (but  the  plane  stress 
case  can  be  derived  by  a  trivial  change  in  the  values  of  the  elastic 
constants).  Two  different  kinds  of  boundary  value  problem  are 
considered. 

In  Problem  (a)  (treated  in  S 2 ,1 )  we  assume  that  the  surfaces  of 
the  strip  are  constrained  in  Buch  a  way  that  the  normal  component  of  the 
surfaoe  displacement  and  the  surfaoe  shearing  stress  both  vanish;  this 
problem  is  exactly  that  of  determining  the  distribution  of  stress  in  an 
infinite  two-dimensional  elastic  medium  containing  an  infinite  row  of 
identical  Griffith  cracks  equally  spaced.  This  boundaiy  value  problem  is 
solved  by  reducing  it  to  the  solution  of  a  pair  of  dual  series  relations 
which  is  in  turn  shown  to  be  equivalent  to  that  of  an  integral  equation 
which  can  be  solved  easily  for  an  arbitrary  distribution  of  internal 
pressure.  In  §2.2  the  solution  in  the  case  where  the  internal 
pressure  is  constant  is  considered  in  more  detail;  the  critical  value  of 
the  internal  pressure  is  calculated  using  Griffith’s  criterion 
^formulae  (2.35)  and  (2.3 6)  belowj  and  the  effeot  of  the  fact  that  the 
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strip  is  of  finite  width  is  illustrated  by  calculating  the  variation 
with  a/c  of  the  percentage  ohange  in  the  value  of  the  critical 
pressure  from  the  value  in  the  'infinite'  case  c  »  a  (Table  1  and 
Pig.  2  below).  The  shape  of  the  crack  and  the  variation  of  the  normal 
component  of  stress  across  the  line  of  the  crack  are  also  shown  for  some 
non-zero  values  of  the  ratio  a/c  (Pigs.  4  and  5  below). 

In  Problem  (b)  (considered  in  §2.2)  it  is  assumed  that  the 
surfaces  of  the  strip  are  free  from  applied  stress.  Here  the 
analysis  is  more  complicated}  the  problem  is  again  reduced  to  that  of 
solving  an  integral  equation  -  in  this  instance  a  Predholm  equation  of 
the  second  kind  -  but  the  equation  does  not  appear  to  have  an  exact 
(analytical)  solution  so  that  in  any  given  case  it  would  have  to  be 
solved  numerically. 

The  discussion  in  §3  is  an  immediate  generalization  to  the  case 
of  axial  symmetry  of  the  methods  used  in  §2  so  that  the  problem  is  now 
that  of  determining  the  distribution  of  stress  in  a  very  long  circular 
cylinder  which  has  an  internal  penny-shaped  crack  which  has  its  centre 
lying  on  the  axis  of  the  cylinder  and  its  plane  normal  to  that  axis. 

The  stress  field  in  the  cylinder  is  due  to  the  application  of  pressure 
to  the  surfaces  of  the  crack.  Again  two  kinds  of  problem  axe 
considered. 

In  Problem  (a)  (discussed  in  §3.1  and  by  an  alternative  method  in 
§3  .2)  it  is  supposed  that  the  shearing  stress  on  the  surface  of  the 
cylinder  and  that  the  radial  component  of  the  surface  displacement  are 
both  identically  zero.  For  an  arbitrary  (but  axisymmetric)  distribution 
of  pressure  on  the  crack  surface  the  determination  of  the  various 
quantities  of  physical  interest  is  made  to  depend  on  the  solution  of  a 
Predholm  integral  of  the  second  kind.  An  iterative  solution  of  this 
equation  with  the  free  term  corresponding  to  a  constant  internal 
pressure  and  valid  for  small  values  of  the  ratio  a/c  is  given  in  §3.3; 
for  values  of  a/c  lying  between  1  and  2  it  is  necessary  to  solve  the 
question  numerically  and  numerical  values  of  the  unknown  function 
corresponding  to  a  set  of  values  of  a/c  are  reported  in  the  same 
section.  Prom  this  solution  it  is  a  simple  matter  to  calculate  the 
variation  with  a/c  of  the  critical  pressure  required  to  cause  the 


3 


crack  to  spread  and  that  of  the  stress  intensity  factor  introduced  by 
Irwin;  these  are  shown  in  Pigs.  9  and  10  respectively.  The  results 
so  obtained  are  compared  (in  Pigs.  11  and  12)  with  those  obtained 
in  the  analogous  case  in  plane  strain  (§2.1  above)  and  it  is  shown 
that  as  far  as  an  "engineering"  approximation  is  concerned  the 
size  effect  in  the  axisymmetric  case  can  be  simply  gauged  from  that  in 
the  plane  strain  case. 

Finally  in  Problem  (b)  ( treated  in  §3.4-)  it  is  assumed  that  the 
curved  surface  of  the  cylinder  is  free  from  applied  stress. 

Again  the  problem  of  determining  the  quantities  of  physical  interest 
corresponding  to  an  arbitrary  axisymmetric  distribution  of  stress 
on  the  crack  surfaces  is  reduced  to  that  of  solving  a  Fredholm 
integral  equation  of  the  second  kind  but  now  the  actual 
calculations  are  more  complicated  since  the  kernel  of  the  integral 
equation  is  a  function  of  the  Poisson's  ratio  of  the  material 
of  the  cylinder.  Again  an  iterative  solution  and  numerical 
solutions  for  a  constant  internal  pressure  (and  for  Poisson’s 
ratio  equal  to  i)  are  given  and  the  variation  with  a/c  of  the 
critical  pressure  and  the  Irwin  stress  intensity  factor 
calculated  (Of.  Pigs.  l6,  17).  These  two  diagrams  also  afford 
a  comparison  between  the  solutions  of  Problems  (a)  and  (b);  if 
a  c  ^  c  there  is  little  difference  between  the  size  effects 
in  the  two  cases  but  if  a  >  £  c  the  effect  is  more  pronounced 
in  the  case  in  which  the  radial  component  of  the  surface  dis¬ 
placement  of  the  cylinder  is  zero. 
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2.  Distribution  of  Stress  in  a  Strip  of  Finite  Width  containing  a 
Griffith  oraak. 


We  begin  by  considering  the  distribution  of  stress  in  a  very 
long  strip  of  uniform  width  2  o  which  has  a  Griffith  crack  of  length 
2a  in  its  interior.  We  shall  assume  that 
the  oraak  is  perpendicular  to  the  edges  of 
the  strip  and  that  its  centre  lies  on  the 
central  line  of  the  strip  (Of.  Fig.  l)  and 
that  the  state  of  stress  in  the  strip  is 
due  to  the  application  of  pressure  of 
prescribed  value  to  the  surfaces  of  the 
crack.  In  addition  we  shall  make  the  usual 
assumptions  of  the  classical  (infinitesimal) 
theory  of  elasticity  and,  in  particular,  we 
shall  solve  the  equations  corresponding  to 
a  state  of  plane  strain  in  the  strip.  The 
results  in  the  case  of  plane  stress  can  easily 
be  deduced  by  a  trivial  change  in  the  values 
of  the  elastic  constants.  We  employ  the 
notation  of  Green  and  Zeraa  ( 1 954)  . 

We  shall  consider  two  types  of  boundary  value  problem.  In 
problem  (a)  we  assume  that  the  surfaces  of  the  strip  are  constrained 
in  such  a  way  that  the  normal  component  of  the  surface  displacement 
and  the  shearing  stress  both  vanish.  This  problem  has  been 
considered  by  Westergaard  (l  93  )  and  Green  and  England  (l  963)  -  in 


the  form  of  a  discussion  of  the  state  of  stress  in  an  infinite  thin 
plate  containing  an  infinite  row  of  identical  Griffith  cracks  evenly 
spaced;  the  solution  given  here  is  that  derived  recently  by  Sneddon 
and  Sriyastav  (1963)  whioh  gives  the  results  in  a  form  suitable  for 
numerical  calculation  and  has  the  added  advantage  that  it  can  be 
generalized  to  provide  the  solution  of  the  analogous  problem  in 
three  dimensions  (Of.  §3.2 below).  In  problem  (b)  we  assume  that 
the  surfaces  of  the  strip  are  free  from  stress.  This  problem  does 
not  seem  to  have  been  discussed  previously;  again  we  follow  the 
method  reoently  devised  by  Sneddon  and  Srivastav  (1963)  which  is  a 
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generalization  to  the  elastost&tio  case  of  a  method  used  "by  Sneddon  (1962) 
in  the  solution  of  a  simple  problem  in  electrostatics. 


2.1 .  Solution  of  Problem  (a) . 

In  the  discussion  of  the  first  of  these  two  problems  we  take  our 
unit  of  length  to  be  such  that  the  width  of  the  strip  is  2  jr  (i.e.,  in 
the  notation  of  Pig.  1  we  take  c  =  *r  ).  The  problem  of  determining 
the  stress  in  the  neighbourhood  of  the  crack  is  then  obviously 
equivalent  to  that  of  determining  the  stress  in  the  semi-infinite 
strip  -i  s  y  «  1  ,  x)  0  when  the  boundary  x  =  0  is  subjected  to  the 
conditions 


=  0, 

-  V  i  y  «  7T  , 

(2.1) 

=  -2  \x  f  (y)  , 

-  a  <  y  <  a, 

(2.2) 

=  0, 

a  <  |  y  |  «  ir 

(2.3) 

where  jj  denotes  Lame's  constant,  and  the  function  f(y)  is  prescribed* 
If  we  assume  that  the  edges  of  the  strip  are  constrained  in  such  a  way 
that  the  normal  component  of  the  displacement  and  the  shearing  stress 
both  vanish  then  we  must  ensure  that  when  y  =  +  n  , 


a 

xu 


=  0, 


0  <  x  <  00  . 


(2.4) 


We  further  assume  that  as  x  — ^  »  the  components  of  stress  and 
displacement  all  tend  to  zero.  We  shall  satisfy  this  condition  by 
finding  a  solution  which  satisfies  the  conditions 


(7 

xu 


u 

X 


=  0, 


-ir  5  y  <  7r 


(2.5) 


on  the  line  x  =  6 


and  then  finding  its  limiting  form  as  6  — ^  » 
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If  we  assume  expressions  of  the  form 


CO 

\  uq  ^  ^  aR  ooseoh(n  5)  fn(  6  -  x)cosh  n(  6  -  x)  - 

'  '/  n»o 


-  (2  -  2  n  +  n5  coth  n  5  )sinh  n(  6  -  x)J  oos(ny) 


(2.6) 


QO 

u  =  -  )  a  cosech(n  6)1  n(6  -  x)sinh  n(  6  -  x)  +  (l  —  2  rj  -  n  6cothn6)  x 

y  Z_/  n  L 

n=  i 

x  cosh  n(  6  -  x)J  sin(ny)  (2.7) 


for  the  components  of  the  displacement  then  (if  we  are  considering  plane 
strain)  the  components  of  stress  are  given  by  the  equations 


axx  0  “  n  )u 

—  mm  m 


-  T)  )u  c->  r 

...I ...  _  \  n  a  coseoh(n  6  )  (l  +  n  6  coth  n  6  )  x 

2m  2(l  -  2v  )  6  "  L 

x  cosh  n(6  -  x)  -  n(6  -  x)sinh  n(  6  -  x)  J  cos  ny 


guu  "u. 


2g  2(1  -  2  r?) 6 


OO 

-  n  a^  cosech(n  6)  ^(l  -  n  6 coth  n  6)  x 


n»i 


oosh  n(  6  -  x)  +  n(  6  -  x)sinh  n(  6  -  x)j  oos(ny) 


2m 


a  ^  n*  an  °osech(n  6  )  ^ 


oosech(n  6  )  (  6  -  x)cosh  n(  6  -  x)  -  6  coth  n  6  x 


n»i 


sinhn(  6  -  x)J  sin(ny) 
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where  77  denotes  Poisson's  ratio. 

This  solution  obviously  satisfies  the  condition  (2.1 )  on  x  =  0, 
the  conditions  (2.4)  on  y  =  +  ir  ,  and  the  conditions  (2.5)  on  x  =  6  . 
Also  on  x  a  0  we  have  the  expressions 


00 

u  =  2(l  -  77  )  fw  a  +  )  a  cos  nyl 

*  l  °  u  "  J 

n*  1 

00 

=  -  %  a  &o  -  ^  nan|^1+  k(n6)J  00s  ny 


where 


u. 


a  = 


0  2(1  -  n  ) 


a  = 


2(l  -  7)  )2 
(1  -  2  7)  )6 


(2.3) 


and 


k(5) 


g+  2e"€  sinh£ 
sinh*  g 


(2.9) 


Hence  we  shall  have  solved  the  relevant  boundary  value  problem  if  we  can 
find  a  sequence  of  constants  J  an  {  satisfying  the  dual  series  relations 


a0  + 


n  an  |l  +  k(n  6  )j  cos(ny)  =  f(y),  0  ^  y  <  a,  (2.1 


0) 


<0 

^a0  +  a n  cos(ny)  =  0,  a  <  y  5  w  .  (2.11 ) 


n*i 


If  to  now  consider  the  case  of  a  long  strip  so  that  b  »  ir  we 

find  that  equations  (2,6)  and  (2.7)  reduce  to 

m 

us  ■  &  u0  6  m*)  +  [..  2(1  -  n  )  +  nxJe‘**oos(ny), 

'  n«i 


(2.12) 
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00 

ug  =  -  an  j(l  -  2  r?)  -  nxje“nx  sin(ny), 
n»  i 


(2.13) 


and  that  the  corresponding  components  of  stress  are  given  by  the 
equations 


2m 


f  »a  - 


00 

^  n  a  n  (l  +  nx)  e“  ns  cos(ny), 


n*  i 


(2.14) 


2m 


oo 

£■ 


a  n  (l  -  nx)  e“  1  cos(ny). 


n*  i 


(2.15) 


fjUL  a  . 

2m  ' 


00 

I" 


2  a  e”"1  ain(ny), 


n«  i 


(2.16) 


where  the  constant  /3  is  defined  by  the  equation 


p  =  ~-ul 

(l  -  2  77  )  6 


and  the  constants  jar  satisfy  the  dual  series  relations 


00 

£  «  a  o  +  l  n  an  cos(ny)  =  f (y)  ,  0  «  y  <  a, 


n  =  i 


(2.17) 


i  *0  +  ^  an  oos(ny)  =  0, 

n*i 


a  <  y  «  it. 


(2.18) 


In  the  case  of  an  infinitely  long  strip  we  may  take  a  =  0. 

To  solve  the  pair  of  dual  series  equations  (2.17)  and  (2.18)  we 
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make  use  of  a  method  due  to  Srivastav  (l963)»  We  make  the 
assumption  that  when  0  <  y  <  a 


_ bOO.^-- 

d/( cos  y  -  cos  t) 


(2.19) 


and  reduce  the  problem  to  that  of  determining  the  function  g(t).  Prom 
the  theory  of  Fourier  series  it  follows  that 


a 


0 


g(t)  dt, 


I  g(t)|Pn(cost)  +  Pn_i(cos^| 


dt 

(2.20) 


where  Pn  denotes  the  Legendre  polynomial  of  degree  n. 

Now  if  a  =  0  In  equation  (2.17)  we  may  integrate  both  sides  of 
the  equation  to  obtain  the  relation 


^  a  n  sin  ny  =  P(y), 
n*i 


0  «  y  <  a 


(2.21) 


where  the  function  P(y)  is  defined  by  the  equation 

F(y)  =  f  f(u)  d  u.  (2.22) 

*  0 

If  we  substitute  from  equations  (2.20)  into  the  left  hand  side  of 
equation  (2.21 )  and  interchange  the  order  of  integration  and  summation 
we  find  that  equation  ( 2 . 2l )  is  equivalent  to  the  integral  equation 


[  S(y,  t)g(t)  d  t  =  F(y),  o  <  y  <  a  (2.23) 

*  o 

where  the  kernel  S(y,  t)  is  defined  by  the  relation 

* 

8(y»  t) 

[p"<o°**>+  Pn-i(oo«tJein(ny)  (2.24) 

riM  V 
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It  is  easily  shown  that  the  function 

oos(£  y)H(y  -  t) 

t)  =  -  (2.25) 

/t/(cos  t  -  cos  y ) 

has  as  its  Fourier  half -range  sine  series  in  the  range  0  $  y  <  jr  the 
expansion  on  the  left  hand  side  of  equation  (2.24) .  Substituting  this 
expression  for  the  kernel  into  the  integral  equation  (2.23)  we  find 
that  it  reduces  tc  the  form 

u  g(t)  oos(£  y)  dt 

-  F(y),  0  $  y  <  a  <  v  (2.2 6) 

0  a/ (coa  t  -  cosy) 


which  can  easily  be  shown  to  possess  the  solution 

2d  m  F(y)sin(£  y)  dy 


g(t) 


*  a  rx 
ir  dt  J  o 


/V/(cos  y  -  cost) 


(2.27) 


where  F(y)  is  given  by  equation  (2.22). 

From  equation  (2.12)  we  have  that  the  surface  displacement  of 
the  crack  u^  (y)  =  u^  (0,  y)  is  given  by  the  equation 

00 

(y)  =  2(l  -  n  )  ao  +  ^an  cos(ny)  J  ,  0«  y 


<  a 


n*i 


so  that  it  follows  from  equation  (2.19)  that 

u(t0)  (y)  =  2(1  -  n  )  cos(£  y)  f  - - *  t 

<»  y  fv(ooB  y  -  cos  t) 


(2.28) 
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Similarly  the  normal  component  of  stress  across  the  line  of  the 
crack  (y)  =  a  sJt  (0,  y)  is  given  by  the  equation 


ffxx  (y)  =  -2n 


E- 


an  cos(ny)  =  -2 tu 


JL  j  ± 

9y  t2 


a  aim 


n*  i 


and  substituting  from  equation  (2.19)  we  find  that  this  is  equivalent 
to  the  equation 

=  -2u  fcos(£  y)  f  - siPL^t -  1  #  (2.29) 

9y  v,  J  0  V(cos  y  -  cos  t)  J 


The  other  physical  quantity  in  which  we  are  interested  is  the 
strain  energy  of  the  crack  in  the  case  in  which  the  normal  stress 
across  the  surface  of  the  crack  is  constant  .  This  is  easily 

shown  to  be  given  by  the  equation 

W  =  -2  (‘  »«(y)uf(y)dy 

0 


If  we  substitute  the  expression  (2.28)  for  u  °  (y)  into  this  equation 

and  interchange  the  order  of  the  integrations  we  find  that 


4A/2(l  -  h2)pp 

E 


cos(^  u )  d  u 

— - (2*30) 

4^(cos  u  -  cos  t) 


2«2.  Crack  opened  out  by  constant  internal  pressure. 

If  the  crack  is  opened  out  by  a  constant  internal  pressure  p 
then,  in  the  notation  of  equation  (2.2)  we  take  f(y)  =  p  /2  /i .  It 


> 
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follows  from  equation  (2,27)  that  in  this  ease  we  need  to  evaluate 
dl/dt  where 

x  sin(£  *)  dx 


■/. 


o  ^(cos  X  -  cos  t) 


or,  what  is  the  same  thing, 


ft  =  -  £  ain^  t)  ft  » 


dl 


dt 


dr 


where  r  =  cos  ^  t  and 


_  2  A/2P  IsgsTI-zkg,. 

Jt  //(z2  -  r2) 


Using  the  formula  for  integrating  by  parts  we  find  that 
log(z  +  V(z2  -  T  Z)]dz 


'  2  V2{  /r 


V(1  -  ■  *) 


=  cos " 1 r 


,  log  rj 


from  which  it  follows  that 


ft  =  “  =  "  V2  ir  sec(£  t) 


and  henoe  that 


ft  “  $2  *)• 


From  equation  (2.27)  we  can  therefore  deduce  that 

Po 


s(t)  = 


tan(£  t) 


(2.31 ) 


The  surface  displacement  (y)  of  the  cr&ok  is  then  given  by 


equation  (2.28)  in  the  form 

2  4/2(l  -  r?2)p  f 

u(  ^  (y)  = - -  cos(£  y)  J 


E 


a  tan(4r  t)  d  t 
y  V(cos  y  -  00a  t) 


The  integration  is  elementary  and  it  is  readily  shown  that 
Co)  '  '  4(l  -  r)*)pn  fcos(£  y)  +  V(°os24-y  -  cos*  •§•  a) 


(y)  = 


E 


log 


l- 


cos  g  a 


] 


|y!  < 


(2.32) 


Similarly  the  normal  component  of  stress  across  the  line  of  the  crack 
is  given  by  equation  (2.29)  in  the  fora 


a ^**(y)  =  “2po  —  j^cos  ^y  J 


a  -g  sec(|-t)  tan(^t)  dt 
o  >/(  sec  2  -g-  y  -  sec  2  ^  t 


] 


a  <  |  y  |  in 


,2po 


[sln- 


(sec£  t  cos£  y) 


t=a 

t=o 


=  Pn 


' _ 

*.//  * 


ain(£  y) 


'/(cob  £  a  -  oos2  |r  y) 


-  l]  ,  a  <  |y|  <  jr  .  (2.33) 


W  = 


In  a  similar  way  equation  (2.30)  yields  the  formula 
4  V2(l  -  >?2)pf  /•  a  ft  oos(^y)  dy 


E 


j  tan(|-  t)dt  f 
•>  0  *  < 


0  //(cos  y  -  oos  t) 


14 


for  the  el&stio  energy  of  the  oraok.  Now 
.  t  cos(£  y)  dy  1  ,  t  cos(£  y)  dy  v 

J  o  ^(cos  y  -  cos  t)  H/2  I  o  ^(sin2  £  l  —  sin2  •j  y)  ~  V2 


so  that  in  our  system  of  units 


W  = 


4(l  -  ha)  ffp2  f  a  8(l  -  h2)p2tr 


i  /a  oil  -  n  ;p  ir 

■*  I  tanC’j  t)dt  =  ■■■  — .  ?  ■  log(sec^-a)  . 

Jo  E 


In  conventional  units,  if  2o  is  the  width  of  the  strip 

!<r  z  )  • 


8(1  -  rj*)p2  o* 

W  =  .  log/ 


irE 


(2.34) 


If  the  surface  tension  of  the  material  is  T,  the  surface  energy  of  the 
crack  is  U  =  4a  T  and  Griffith's  criterion 


£  (*-H)=0 

for  the  critical  value  pcr  of  p^  when  the  crack  length  is  2  a  leads 
to  the  equation 


p  .  f  — S —  cot  as  ] 

cr  (l  -  r? 2  )o  2c 


In  the  infinite  oase  (i.e.  o  = »  }  the  corresponding  value  is  known  to  be 

£ 


O)  _  f  2  BT  'j 
°r  =  l  »(1  -  r,‘)J 


so  that  we  may  write  the  last  result  in  the  fprm 

(«•) 


Per  =  Pcf  w 


(2.35) 


where 


a)  (z)  =  ^irz  cot(£irz)j 


(2.36) 


If  we  denote  by  n  the  percentage  change  in  the  value  of  the 
critical  pressure  pQr  from  the  value  p^”^  corresponding  to  the 
infinite  case  (o  »  a)  so  that 

C00)  _ 

n(a/c)=  — — £se  x  i  oo 


n  (a/o)  =  100  Jl  -  cot  gj  2  J 


(2.37) 


The  variation  of  the  functions  ui  and  II  with  a/c  is  shown  in 
Table  1  and  that  of  n  is  shown  graphically  in  Pig.  2. 

We  now  return  to  the  consideration  of  the  shape  of  the  crack  and 
the  value  of  the  normal  component  of  stress  across  the  line  of  the 
crack.  Putting  y  equal  to  zero  in  equation  (2.32)  we  find  that 
u^  (o)  =  (  ,  the  depth  of  the  crack  is  given  by  the  equation 


e  =  d(s/c) 


(2.38) 


where  e'  '  denotes  the  depth  of  the  crack  in  the  case  c  »  a  and  is 
given  by  the  equation 


,—,  2(1  -n‘)p0. 


(2.39) 


(in  conventional  units)  and  the  function  d(z)  is  defined  by  the  relation 


d(z)  =  log  tan  j  £  ir  (z  +  1 )  j  . 
V  Z  V.  J 


(2.40) 


I 


3.1 

3.2 

3.3 

3.4 
3.3 

3.6 

3.7 
0.8 
0.9 
0.92 
0.94 
0.96 
0.98 
0.99 


0) 

n 

0.9958 

0.42 

0.9833 

1.67 

0.9617 

3.83 

0.9300 

7.00 

0.8862 

11.37 

0.8274 

17.26 

0.7489 

25.11 

0.6389 

36.11 

0.4732 

52.68 

0.4272 

57.28 

0.3735 

62.65 

0.3080 

69.20 

0.2256 

77.44 

0.1562 

84.38 

Figure  2 


Fig.  2  The  variation  with  a/c  of  rr,  the  percentage  increase  in  the 
critical  pressure  over  its  value  in  the  case  c  =  °°,  for  constant  internal 
pressure. 
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Values  of  d(z)  for  a  range  of  values  of  z  between  0  and  1  are  listed 
in  Table  2,  and  the  variation  with  s/e  of  e  ,  the  depth  of  the  era  ok, 
for  a  fixed  internal  pressure  pfl  is  shown  in  Pig.  3. 

Written  in  terms  of  conventional  units  equation  (2.32)  takes  the 

form 


X°) 


(y)  =  « <wi> .  ff :  ^8  (o°s(  §: )  +  V(°ob*  g  -  cos2g  jj 


-  log  cos 


<8>j- 


(2.41) 


The  variation  of  the  ratio  u^  (y)/e^  with  y  for  four  values  of  the 
ratio  a/o  is  shown  in  Pig.  4;  the  curve  corresponding  to  the  value 
a/c  ~  0  is  an  ellipse 

u«  (y)  =  (1  -  y2/a2)* 


and  it  can  be  seen  from  the  curves  of  Pig.  4  that  even  a  substantial 
increase  in  the  ratio  a/c  does  not  appreciably  affect  the  shape  of 
the  curve  although,  (as  we  should  expect  on  physical  grounds)  the 
’minor  axis'  increases  by  nearly  8C0S  as  a/c  increases  from  0  to  0.9. 

Similarly,  it  follows  from  equation  (2,33)  that  (in  conventional 
units)  the  normal  component  of  stress  aoross  the  line  of  the  crack  is 
given  by  the  equation 


*«(y)  =  p0  [8in(  §  )  [°os2  (  8  )  _  C08‘(  S  [I  ^_1]  » 

a  <  |y|  <  c.  (2.42) 


The  form  of  the  variation  of  a  ^°'(y)  with  y  in  three  cases  is  shown 

xx 

in  Fig,  5,  (the  same  value  of  po  being  ohosen  in  each  case).  As  we 
should  expect,  we  find  that  the  shape  of  the  curve  is  muoh  the  same 
whatever  the  value  of  the  ratio  a/ c  but  the  intensity  is  greater  the 
higher  the  value  of  s/o. 


Figure  3 

Fig.  3  The  variation  with  a/c  of  the  depth,  e,  of  the  crack  in  the 

( QQ  \  2  , 

case  of  constant  internal  pressure  pQj  e  *  2(l-r\  )pQa/E  is  the. depth 
in  the  case  a/c  ■  0. 


Fig*  5 

across 


(0)  (v)  the  normal  component  ol  stree 
The  variation  «ith  y  •'% J  constant  internal  preaanre  P( 

the  line  o£  the  orach  in  the  case 
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The  stress  intensity  factor  N  is  defined  by  the  equation 


N  s  lim  V?  a  0  (a  +  r) 

XI 

r  o  + 


=  p. 


If  o  >  >  a  N  takes  the  value 


=  P„(*«)i 


so  that  we  may  write  the  percentage  increase  of  N  over  N 

M 

- 7ST~  =  n(a/c) 

IP  ; 


6») 


as 


where 


n(z)  =  [[^  «)J 


-  1 


J 


x  100. 


The  function  n(z)  is  tabulated  in  the  second  column  of  Table  2: 
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2.3.  Solution  of  Problem  (b). 

In  the  second  problem  the  boundary  conditions  (2.l)  -  (2.3) 
remain  but  the  conditions  (2,4)  are  replaced  by  the  conditions 


axu  = 

ffuw 

=  0, 

If  we  denote  by  j* 

s 

transform  operators 

and 

> 

y.j 

f  * ; 

3  J 

1  <l>  i 

w 

y  =  +  c,  0<x<oo 
the  Fourier  sine  and  cosine 


<Kx,y)sin(£  x)dx 


<f>(x,y)coa(&x)d.x 


(2.43) 


then  we  can  write  the  boundary  conditions  (2.43)  in  the  equivalent 
forms 


Uv L 


]  [  C<'»«'Dy=0  **->«]  =  o  c^) 


We  reduce  the  solution  of  this  problem  to  that  of  a  Fredholm  integral 
equation  of  the  second  kind  by  a  method  similar  to  that  employed  recently 
by  Sneddon  ( 1 962)  in  the  solution  of  a  boundary  value  problem  of  mixed 
type  in  eleotrostatios.  We  assume  a  stress  field  of  the  form 


u 

x 

X 


=  /  [[fte)  +  2(l  -n)g(5i)°osh(gy)  +  fcygte) 

sink(5y)Js  ’  sin(gx)d5  +  ^  j  5"1  4>(C)  x 


x 


X  (2(1  -r?  )  +  c  *3  e“  Cx  oos(  5  y)d£  , 


(2.45) 
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Uy  =  '  ■LfVT  /  jjffe)  -  (1  -  2n)g(z}]  sinhfey)  +  £yg(g) 


x  cosh(  g  y) 


5’1  cos(5x)dg  S"V(C)x 


xQ  -  2ri  -  C  x)e~  ^ X  3±n(  £  y)  d  £  , 


(2.46) 


where  E,  as  above,  denotes  Young's  modulus,  n  denotes  Poisson's  ratio 
and  f ,  g,  and  <f>  are  functions  to  he  determined  from  the  boundary 
conditions.  If  we  substitute  these  expressions  into  the  stress-strain 
relations  we  find  that  the  corresponding  stress  components  are 


xx 


=/?/:  j(f  +2g)cosh(£y)  +  £ygsinh(£y)J  cos(£x)d£ 


-JSI  $(0(1  +  S  x)e"  ^  x  oos(£  y)d£  , 
r  TT  J  0 


(2-47) 


a 

uu 


/  —  j  jjf  cosh(gy)  +  £  ygsinh(Cy)  cos(gx)d 


-vT/. 


~  ^  <#>(C)(l  -  S  x)e~^x  oos(C  y)d£  , 
n  J  o 


(2.48) 


xy 


■v?r  £(f  +  g)sinh(gy)  +  gygcosh(Cy)J  sin(gx)dg 

*  ^’x  f  £  *(£  )  e"  C  x  sin (£  y)d£  .  (2.49) 

^  #  0 


This  solution  therefore  has  the  following  properties:- 
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On  x  =  0;  uxS^-Y^'yi J  £y)dg  > 


=</?/;  ^  (f  +  2g)cosh(£y)  +  £yg  sinh(£  y)j  d£  - 

-y^l'dy/'  C1  ^(C)sin(C  y)d£  , 

•  0 


axU  =  °> 


=  ■-  o:  3^  ^cjTyy  J  j  =  -  f  oosh(gc)  -  g  eg  sinh(£c)  - 

.bSir  g.^(g)coa(g.o)de 

Jo  u2  +  c2)2 

?s[axU  »  =  +  ^(f  +  g)sinh(gc)  +  5 cgcosh(  gc)J  - 


-  A£  j  £2  <dc)aintio)dC 

+  n  Jo  {V  +  £2)2 

It  follows  from  these  equations  that  the  boundary  condition  (2.1 ) 
is  automatically  satisfied  and  that  the  remaining  conditions  (2.2)  and 
(2.3)  on  x  =  0  are  satisfied  if  we  choose  f ,  g,  <f>  such  that 


*  U  \ 

_  VI  jL  f 

y  n  dy  Jo 


r* 

C'tiC  )cos(5  y)d£  =  0,  a  <  y  $  c 

j  0 


[  de  - 

0  <  y  <  a. 

(2.50) 

a  <  y  s  c 

(2.51) 
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where  we  have  written  p(y)  for  2  tu  f (y) .  Similarly  the  conditions 

(2.44)  on  y  r  +  c  will  be  satisfied  if 

f  cosh(g  o)  +  £  o  g  sinh  (5  o)  -2  i  (€)  (2.52) 

(f  +  g)  sinh(£  o)  +  5  c  g  oosh  (£  o)  =  2i2(€)  (2.53) 

where  the  functions  i^(£) ,  ig(£)  are  defined  by  the  integrals 

i  <«)  =  ruttiaalSstoi  ,  la(5)c2S  s  (2-54) 

IT  •'o  (€  +  £  )  V  (£2  +  C 2)2 

If  we  make  the  representation 

*(£)=£  fa  t*(t)J0(£t)at  (2.55) 


of  the  function  <f>  ($)  then  it  is  easily  seen  that  the  equation  (2.51 ) 
is  satisfied  whatever  the  form  of  the  function  tp  (t)  but  that 
equation  (2.50)  is  satisfied  only  if 


V(y8- t2) 


-  f  (f  +  2g)cosh(£y)  +  £y  g  sinh(€y)j 

J  0  v.  J 


d€  = 


=  /l  p(y)» 


0  s  y  <  a. 


If  we  integrate  both  sides  of  this  equation  with  respect  to  y  from  0 
to  y(  <  a),  we  see  that  it  is  equivalent  to  the  equation 


/  V(&*t^<t*)  "/  +«)8iQh(£y)  +  ysoosh(5y)J  <*£  = 

=  p(y), 


0  <  y  <  a 
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where 


P(y) 


=  f/\ I  p(t)dt,  0  s  y 


<  a. 


(2.56) 


Regarding  this  as  an  integral  equation  of  Abel  type  for  <l>  (t)  in 
terms  of  known  functions  f ,  g,  P  we  see  that 


2  r«»  dr 

t*(t)  =  -  e-’(f  +  g)as- 

ir  »  o  dt  *  < 


d  r  t  y  sin  h(6  y)*y 


o  y(t2-y2) 


y2cosh(5y)dy  2  d  -tyP(y)dy 
//(t2-y2)  +  Jr  dJ0  ♦'(t2  -  y2)  ’ 


Using  the  results 


2  d  n  y  sinh(€y)  dy 

Tr-etiUt), 


/ 

ir  dtic 


ir  dt  Jo  V(t  -  y  ) 


2  d  rt  y 2  cosh(  £  y)dy  -  .  n 


2  d  /•*  y  P(y)  dy 


ci  Q.  f 

IT  dt  J  < 


0  y(t  -  y 


;  yr/: 


p(y)  dy 

V(t“  -  y 2) 


we  find  that  the  relation  between  ^  ,  f ,  g  and  p  can  be  written  in 
the  form 


<l>  (t)  =/  ^  (f  +  2g)lo(  «t)  ♦get  !,(**)]  dg  + 

/2*  f 1  p(y)  dy 

if  0  </(t*  -  y2) 


(2.57) 


On  the  other  hand  if  we  make  the  substitution  (2  ^>5)  into  the 
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expressions  (2.55)  for  the  integrals  1^(5  )>  i2(S  )  find  that 
these  integrals  can  be  written  in  the  form 


ra  r  a  3i 

i(0=€  /  u|}(u)  i,  du,  i2(0  =  -  u^(u)  — - 

J o  Jo  do 


du 


where 


1  a  f°°c os(£  o)j  (  C  u)  d£ 


1  a  f»cosu  cJJ  I 
i3(g,u,  c)  „ — d  — 

ir  35  Jo  C  + 


S' 


The  integral  occurring  on  the  right  hand  side  can  be  evaluated  by- 
formula  (14)  on  p.45  of  vol.I  of  Krdelyi  (l 954)  to  give 


ij( 5  >  u»  G)  -  2 


1  JL 

dZ 


Ze^  I 


,(€u)j  =  \ 


•  X 


jo  -  Cg  )  Io(g  u)  +  SuI^Su)^ 


so  that 


i1  (5  )  =  |  €  e-0^  J  u^(u)  (l  -  c  g)lQ(g  u)  +  £  u  1^(5  u)  J  du  (2.5fc) 
i2(g)  =  f  €  |  &  u  t/> (u)  £(2  -  c5  )lQ(gu)  +  £u  I^CuJj  du  (2.57) 


Now  if  we  solve  equations  (2.52)  and  (2.53)  we  find  that 

f  (  5  )  =  -  ~  11  *  -  — ^  11  j  jg  0  cosh(  Z  0)  +  sinh(  Z  o)  J  x 

2Zo  +  sinh(2  5o)'*  *- 

i^S)  +  €0  sinh(  Z  o)  i2(5)J 


(2.58) 
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8(5)  - 


2g  o  +  sinh(2£  o) 


jsinh(  gc)!^)  +  cosh(g  o)± £( (2.53) 


If  we  substitute  from  equations  (2.56),  (2.57)  into  equations 
(2.S8),  (2.59)  and  insert  the  resulting  values  of  f(g),  g(g)  into 
equation  (2.50)  we  find  that  the  integral  equation  for  the 
determination  of  the  function  (t)  is 

ifi  (t)  -  /  ^(u)  K(t,  u)  du  =  q(t) ,  0  <  t  <  a  (2.60) 

•  o 


where  the  kernel  K(t,  u)  is  defined  by  the  equation 

K(t,  u)=u  |  - SA-S - f  (2oz5  2  -  6c  5  +  5  +  3e"2og)x 

Jo  2go  +sinh(2£c)  v- 

x  io(«  u)  Io(5t)  +  (3  -  2cg  +  e"2o?)[iu  I#(5t)l  (5  u)  + 

+  5t  I0(5u)  I/St)]  +  2g2utIT(eu)  If(  €t)J  (2.61) 


and  the  free  term  is  defined  by  the  equation 

/T  r  p(y)  ay 

V  w  Jo  y(t2  -y2) 


(2.62) 


In  the  oase  in  which  the  internal  pressure  p(y)  is  a  constant 
p0 ,  say,  we  have  q(t)  =  (-g-  v  p  ,  so  that  if  we  write 


'll  (t)  =  (i»)^P0#,(i/o) 

we  find  that  <p  ^  (t)  is  the  solution  of  the  Fredholm  equation 


(2.63) 
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ra 

0  (t)  -  /  (v)  L(t,  v)  d  v  =  1 ,  0  <  t  <  a 

1  J  o 


(2.6*) 


where  a  =  a/c  and 


L(t,  v)  =  o  K(ot,  cv) 

r  - -  f(25*  -  6c  +  5  +  3e“2c)l  (£v)l  (ct) 

Jo  2C  +  sinh  2  C  L 


=  v 


(3  -25  +  «"2C)[Cv  l0(5t)l(5v)  +  Ct^CCvJl/Ct)^ 

2C*  tv  I, (;▼)!,(  ct)j.  (2.65) 


An  approximate  solution  of  the  integral  equation  (2.66)  is 


0,(t)  =  1  +  f  L(t,  v)dv. 
*  0 


(2.66) 


Using  the  results 


fa  al(Ca)  /■« 

/  vl(5v)dv  = - -  ,  [  v*I1(5v)dy  = 

Jo  £  J  a 


a*  1,(5 <0 

s 


(2.67) 


we  find 


£  I(t,  >)!-  =  .  £  2;v4Lh  - «« ♦  5  -  so-2') x 

X  I  (Ca)  I  (ct)  +  (3  -  25  +  e“2S  )x 
1  0 

x[jolo(5t)la(5a)+5t  I#(5o)^,(Ctf]  +  252tal2(5a)l,(5t)j. 

(2.68) 
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3. 


Distribution  of  Stress  in  a  Cylinder  containing  a  Penny-Shaped  Crack. 


We  shall  now  consider  the  distribution  of  stress  in  an  infinitely 
long  circular  cylinder  which  has  a  penny-shaped  crack  in  its  interior 
with  the  centre  of  the  crack  lying  on  the  axis  of  the  cylinder  and  its 
plane  perpendicular  to  that  axis. 

(Cf.  Pig.  6).  We  shall  suppose  that 
the  deformation  in  the  cylinder  is 
produced  by  the  application  of  the 
pressure  to  the  surfaces  of  the  crack. 

As  in  the  two-dimensional  case  we  shall 
consider  two  types  of  boundary  conditions 
on  the  curved  surface  of  the  cylinder. 

In  the  first  problem  -  denoted  by 
problem  (a)  here  -  we  suppose  that  the 
shearing  stress  on  this  surface  and  the 
radial  component  of  the  surface  dis¬ 
placement  both  vanish;  this  problem 
has  recently  been  discussed  by  Sneddon 
and  Tait  (1963)  and  we  shall  follow  that 
solution  here.  In  the  second  problem  - 
called  problem  (b)  below  -  we  assume  that 

the  curved  surface  of  the  cylinder  is  free  from  stress;  this  problem  has 
been  considered  recently  by  Collins  (1962)  and  by  Sneddon  and  Welch  (1963). 
We  shall  outline  the  method  of  the  latter  paper  here. 

3.1.  Solution  of  Problem  (a). 

If  we  take  the  radius  of  the  crack  to  be  our  unit  of  length,  and  the 
radius  of  the  cylinder  to  be  o  (  >  l),  then  the  problem  of  determining  the 
distribution  of  stress  in  the  neighbourhood  of  the  crack  is  equivalent  to 
that  of  finding  the  distribution  of  stress  in  the  semi -infinite  cylinder 
p  <  c,  z  *  0  when  its  plane  boundary  z  =  0  is  subjected  to  the  conditions 


Fug.  6 


— 

0 

II 

N 

tlQ' 

0  3  p  <  0, 

(3.1) 

azz  =  ~2^f(p)»  <  U 

(3.2) 

nz  *  °» 

1  <  p  $  0 

(3.3) 
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where  u  ia  Lame's  constant  (rigidity  modulus)  and  the  function  f(p)  is 
prescribed.  If  we  assume  that  the  cylindrical  surface  is  kept  fixed  then 
on  p  =  c  we  must  have 

Up=  a pz  =  0,  0  ^  z  <  oo  (3.4) 

It  is  well-known  03f.  Sneddon  (l95l)>  p.505^  that  a  solution  of 
the  equations  of  elastic  equilibrium  in  the  axially  symmetric  equation  is 
given  by  the  equations 


1  3  r 

u  _  J - fi. 

2  p  9p3z 


(3.5) 


where  x(  P  t  z)  is  an  axisymmetric  biharmonic  function  and  rj  is  Poisson’s 
ratio*  The  components  of  the  stress  tensor  can  be  determined  from  the 
stress -strain  relations;  we  have  • 


A  suitable  type  of  biharmonic  function  for  a  problem  of  this  type  is 
defined  by  the  equation 

X  =  -2/ijV^jAte)  +4(  1  -/])B(s)]l0(£p)  -  epBte)^  (&>)]  sin(?z)dg 

-2p  fV3F(€)(27)  +  5z)e‘&  JQ(5p)d€,  (3.6) 

*  0 


where  A,  B  and  P  are  functions  of  5  alone.  A  solution  of  this  form 
automatically  satisfies  equation  (l.l).  The  corresponding  expressions  for 
0  2Z  and  u  on  the  plane  z  =  0  are  given  respectively  by  the  equations 


zz 


u, 


z 


=  -2 nf  ^(Ate>  "  2nB(^I0(€p)  -5pB(€)l  (€p)jd5 

-2p  f  F(5)J0(5p)d5, 

•  0  w 

=  2(1-0/  r1  P(5)J  (5p)d5. 

Jo  0 


(3.7) 

(3.8) 
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Prom  these  last  two  equations  it  follows  immediately  that  the 
boundary  conditions  (l  .2)  and  (l  .3)  are  satisfied  if  A(fc) ,  B({* ) ,  P(g)  satisfy 
the  dual  integral  equations 


^F(OJ0(5p)d£  +/o  s{[a(0  -  2r?B(S)]l0(5P)  -£  PlKS)^  (5p)  jdg  =  f(p),  0  <p<1, 

[  rM€)J0(£p)<l5  =  0, 

J  0 


1  <  p  ?  c. 


It  is  known  £  Of.  Sneddon  (l960)J  that  the  second  of  these  equations  is 
automatically  satisfied  if  P(^)  is  written  in  terns  of  an  unknown  function 
g(t)  through  the  equations 


P(S)  =5  l  g(t)sin(gt)dt,  g(o)  =0  (3.9) 

*  0 

and  that  if  we  substitute  this  form  into  the  first  equation  of  the  pair 
it  reduces  it  to  the  relation 


rP  g* (t)dt 
0  //(P2  -  t2) 


f(p)  *  A(S)  -  2r,B(£)3l,(6p)  -«pB(Ol(5p)J 


which  can  be  thought  of  as  an  Abel  type  integral  equation  with  solution 

,  s  .  .  2  r  r  .  .  pio(5p)^p  r 

g(t)  =  h(t)  - -/  e  (A-  277B)/  - —-€B/ 

jt^oL  Jo  V(t  -  p  )  0 


•t  p*I  (gp)dp 
-  -  ]  d.  i 

0  v( tz  -  pz)3 


where  h(t)  is  defined  in  terms  of  f(p)  by  the  equation 

2  ft  pf(p)  d  p 


2  ft 

h(t)  =-  / 

w  Jo 


0  //(t2  -  p2) 


(3.10) 


/ 


Making  use  of  the  results 
4  P*  (Sp)dp  sinh(gt) 


0  V(ta-P2) 


ft  p*I  ($p)dp  5t  cosh(£t)  -  sinh(5t) 

j  _ - - - 


o  V(t2  “  P2) 


we  find  that  the  relation  connecting  g(t)  with  A(g)  and  B(g)  may  be  written 
in  the  fora 


g(t)  =  h(t)-i  /  f[A(5)  +  (1  -2n)Bfe3»toh(5t)  -StBfe)  oosh(et)|dg.(3.1l) 
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Further,  if  we  substitute  from  equation  (3.9)  into  equation  (3*8) 
we  find  that  the  normal  component  of  the  surface  displacement  is  given 
by  the  equation 

g(t)  dt 


»,°  (p)  =2(i  -  n)  f'  -■  , 

*  Jp  V(ta  -  p2) 


o  «  p  « 1 


(3.12) 


Another  expression  of  physical  interest  is  the  energy  W  required  to 
open  out  the  crack.  This  is  given  by  the  equation 

W  =  4irp  j  pf(p)u  0  (p)  dp. 

If  we  substitute  the  expression  (3.12)  for  u°  (p)  and  interchange  the 
order  in  which  we  perform  the  integrations  we  find  that 


W=4jrZAi(l-  n)  a  (c"1) 


(3.13) 


where  the  function  co  is  defined  by  the  equation 
o)(c"1)  =  j  h(t)g(t)dt. 

j  0 


(3.14) 


The  expression  (3.13)  for  the  energy  W  is  derived  on  the  assumption 
that  the  unit  of  length  is  the  radius  of  the  crack.  If  the  length  of 
the  crack  were  a  this  expression  would  be  modified  by  multiplying  the 
right  hand  side  by  a3  ,  i.e.  we  should  have 


W  =  4ir2  n  (l  -  n  )  u  (a/o)a3  . 


(3.15) 


In  particular,  if  a  constant  pressure  p  is  applied  to  the  crack, 
f (p)  =  p0/2p  and  so  we  have  ^ 

h(t)=— .  (3.16) 

Wfi 


If  we  write 


we  find  that 


g(t)  =  -4  <(>( t) 

2  4 


8pa  (i  -  n‘) 

W  =  - a3  «  (a/c) 

■*  1 


(3.17) 


(3.18) 


3$ 


where 


u 


i 


The  Griffith  criterion  that  the  crack  may  spread  is 


(3.19) 


_a_ 

3a 


(w  -  u)  =  0 


(3.20) 


where  U,  the  surface  energy  of  the  crack,  is  given  in  terms  of  the 
surface  tension  T  of  the  material  by  U  =  2  waz  T.  This  leads  to  the 
expression 


or 


=  r  ■■■*■?■*■■  iWc) 

1 2(l  -  n  z)a  J  1 


(3.21) 


for  the  critical  value  of  the  applied  pressure  to  cause  the  crack  to 
spread  when  its  radius  is  a.  In  this  equation  the  function  fi  i  (x)  is 
defined  by  the  equation 

n  [x)  =  ^(x)  +  ix«j(x)J  ~s,  (3.22) 

If  c  >  >  a ,  the  critical  pressure  takes  the  value 

_(»)_  f  1  * 

P°r'  l  2(1  -  r\  *)a  J 

so  that  we  may  write  equation  (3.21 )  in  the  form 


Per  =  *2  n  i(a/c)  * 


(3.23) 


Equation  (3. 11 )  gives  one  relation  connecting  the  unknown  functions 
g(t),  A(  €  )>  B(  g  )  and  the  known  function  h(t).  The  two  remaining 
relations  are  given  by  the  conditions  on  the  curved  surface  p  =  o.  It 
is  easily  shown  that  the  values  of  a ^  and  u^  on  the  surface  p  =  o 
corresponding  to  the  form  (3.6)  for  x(  P  >  z)  are  given  by  the  equations 


N, 


=  -2  p  z 
=0  *  0 


J  £F(5)e"5zJ  (Cc)dJ 
*  0 

2  uzj  ej&(5)+2(l  -hjBCg^l/So)  -5cB(€)Io(€c)]  X 

sin(£z)dg  (3.24) 
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[up]  pm  c  =  -j  S'MO  1  -  2n  -  £z  J(Cc)e'Cz  iC 

+ f  jjXfe)  +4(1  -ecB(5)loUc)jcosUz)dg.  (3.25) 

If  we  take  the  Fourier  sine  transform  of  both  sides  of  equation  (3.24) 
and  make  use  of  the  fact,  which  follows  from  equations  (3.4),  that  the 


Fourier  sine  transform  ofja  p 
obtain  the  equation 


p  =  0 


is  identically  equal  to  zero  we 


Afe^Uc)  -  B(£)[c£I0(cg)  -  2(1  -  r?)  1,(05  )]  =  i2  -  i,  ,  (3.26) 


where  the  functions  i  ($),  i  (0  are  defined  by  the  integrals 


,  .  4  p(c)J(cc)  ^  45-  rnc;j  l 

+  C  Tt  J  Q  {£  + 


45*  /*“F(c)J  Uc)ac 


7T 


2V2 


n 


(3.27) 


Similarly  if  we  take  the  Fourier  cosine  transform  of  both  sides  of 
equation  (3.25)  and  make  use  of  the  fact  that  [u0l  vanishes  for 

r J  p  :  C 

all  values  of  z,  we  obtain  the  equation 


A(e)l,tec)  -  B(5)  [  £I0(c  5)  -  4(1  -  n  Jl^c  «)]  =  it  -  ni, .  (3.28) 


Solving  the  equations  (3.26)  and  (3.28)  we  find  that 

--c£in(ce)  ^  io(0  .  .  i.(e) 


A(€) 


1,(0  fc€I0(c€)  ,  ^  1,(0  ,  s  1, 

=  — -  | - 2 - (2  -  „)j  +~L—  ,  B(0  =  — 1 


I  (cO  Ul(c  5) 


I  (cO 


21.(05) 


(3.29) 


Substituting  from  equation  (3.9)  into  the  first  of  equations  (3.27) 
we  find  that  i  (5)  is  related  to  g(t)  through  the  equation 


1(0  =  “  [  g(t)dt  f 
IT  *  0  J  ( 


«  gsin(gt)J  (  gc)d; 
_  2  ,2 

5  +  S 


(3.30) 


The  inner  integral  can  be  derived  easily  from  entry  (5)  on  p.10  of  Vol.2 
of  Erdelyi  (1954)  and  we  find  that 


1,(0  =  ^  lC(o  5)/*  g(u)sinh(5u)du. 

IT  <0 


(3.31 ) 
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Similarly,  it  is  easily  shown  that 

i  (€)  =  “  f  g(u)fsinh(gu)["cEK  (cO+*  (eg)"]  -  £u  cosh(gu)K  (oOdu. 

"  °  1  L  1  1  J  (3-32) 

If  we  substitute  from  equations  (3*31 )>  (3.32)  into  equations  (3.29) 
and  thence  into  equation  (3.11 )  we  find  that  this  last  equation  reduces  to 
the  integral  equation 


g(t)  -  j  K(t,  u)g(u)du  =  h(t) 
J  O 


(3.33) 


in  which  the  kernel  K(t,  u)  is  defined  by  the  equation 
4  f~r 2  I,(o5)K(c5) 


4  rc 

K(t,  u)  =  — 

IT  2  J  O 


IT2  J  O  L  I  (eg) 


sinh  (  gt)sinh  (  gu) 


€*,(<>  g)  r 

+  . .  ■* 1  "  u  cosh(  g u)sinh  (  gt)  +  t  cosh  gt  sinh  gu  1  dg. 

I. (eg)  L 


Now,  using  standard  properties  of  Bessel  functions  (see  Watson  (1944) 
PP. 79-80),  we  have  that 


d 

dg 


K  (eg) 

,  I  (eg) 

u  1 


1  K'(c«)l  (c  g)  -  I* 

=  e  “J - 1 - ; - 1 

J  I*(e  O 


(cg)K  (c  g) 


,  [l0(o€)  +  I2(c  £)]  K  (eg)  +[Kn(cg)  +  Kp(og)]l((cg) 

=  —  -g-  C  .  »  i  'I  i  »■  '  """  1  i 

I'(og) 


1 


€1,(0  5) 

so  that 


sinh  (  g  t)  sinh  (  g  u)  dg 

r«>  d 

X  (eg)” 

L  i*(cg) 

Jo  dg 

I  (eg) 

•  1 

gsinh  (gt)sinh(gu)dg. 


Integrating  by  parts,  we  can  reduce  this  last  term  to 
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/, 


K  (eg) 
0  I  (eg) 


jsinh  (  gt)sinh(  gu)  +  gt  ooah(gt)sinh(  gu)  + 

+  gu  coah(gu)  ainh(gt)"]  dg  • 


We  therefore  find  that  the  kernel  K(t,  u)  of  the  integral 
equation  (3.33)  is  given  by  the  equation 


K(t,  u)=  — 


K  (eg) 

sinh  (  g  t)sinh  (  g  u)  dg 

I  (eg) 


(3.34) 


which  may  be  written  in  the  form 


K(t,  u)  = 


where  the  function  H(  X  )  ia  defined  by  the  equation 
2  r<»  K  (v)  p 

H(  X  )  =  — /  — i—  I  ooah(  X  v)  -  1  dv. 

/Jo  i/v)  L  J 


(3.35) 


(3.36) 


3.2.  Alternative  Solution  of  Problem  (a). 

In  view  of  the  recent  interest  in  the  solution  of  dual  series 
relations  it  is  interesting  to  see  that  problem  (a)  can  be  reduced  to 
that  of  solving  a  pair  of  such  equations  and  that  the  use  of  a  method 
of  solution  due  to  Sneddon  and  Srlvastav  (1963)  leads  to  the  same 
integral  equation. 

If  we  insert  the  biharmonic  function 


*(,,.)  =-2„ 

n*  1 


(2r?  +  Xnz)e"Xnz  Jo(Xnp)  (3.37) 


into  the  equations  (3.5)  we  find  that  it  corresponds  to  an  axisymmetrlo 
displacement  vector  with  components 

n»i 


(3.38) 
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OP 

u,  ■ 2  T, x » %(i  -  ”  * \ z)s'  x"z  ^ 

n*i 


(3.39) 


The  z-components  of  stress  are  then  given  by  the  equations 


)) 


00 

-Z^y  X*  a  (l  +  X  z)e“X"*  J  (X  p), 
/ i  n  n  n  o  " 


n*i 


a  =  -  Zu  z 

3P 


00 

y  X3  a 

U  "  n 

n*  i 


e"  Xnz  J  (Xn  p). 


(3.40) 


(3.41 ) 


It  is  an  immediate  consequence  of  equations  (3.38)  and  (3.41 )  that  the 

boundary  conditions  (3.4)  will  be  satisfied  provided  that  X  ^  ,  X  2,  ...» 

X  ....  are  chosen  to  be  the  positive  zeros  of  the  function  J  ( Xc). 

n  1 

Similarly  it  follows  from  equation  (3.41 )  that  this  choice  of  biharmonic 

function  automatically  satisfies  the  boundary  condition  (3.1 )•  The 

remaining  boundary  conditions  (3.2)  and  (3.3)  are  satisfied  provided  that 

we  can  find  constants  a  to  satisfy  the  dual  series  equations 


OO 


n*  i 


00 

7. 


n*i 


X  aJ(Xp)  =  0, 

n  n  o  n 


1  <  p  «  c 


where  the  X  ^  are  the  positive  zeros  of  J  (Xc),  We  can  readily 
transform  these  equations  to  the  form  considered  by  Sneddon  and 
Srivastav  (1963).  If  we  multiply  both  sides  of  these  equations  by  p 
and  integrate  the  first  with  respect  to  p  from  0  to  p  and  the  second 
from  p  to  c  we  see  that  they  are  equivalent  to  the  pair  of  dual  series 
equatioiis 


\  7 

n»i 


X  n  an  J  (Xnp)  =  F  (p),  0«p<1, 


(3.42) 
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00 

n*  i 


1  <  p  <  o, 


(3.43) 


where  X  ^  ,  X  2 ,  . . . ,  X  n ,  ...  are  the  positive  zeros  of  J  (  X  o)  and 

■P 


F(p)  =  1  f  uf (u)  du. 

J  0 


p  *  0 

It  has  been  shewn  by  Sneddon  and  Srivastav  (1963)  that  the 
solution  of  these  dual  equations  is 

2 


a  = 


n  oV(Xc)  Jo 


j  utt(u)j  (uXn)du, 

*  O 


(3.44) 


(3.45) 


where 


a(p)  =  Y  an  J  (X„p)  =  -  —  I 
l— '  d  p  J  P 


g(t)  dt 


V(t2  -  P2) 


,  0  <  p  <  1  ,  (3.46) 


and  g(t)  is  the  solution  of  the  integral  equation  (3.3)  where  K(t,  u)  is 
the  function  defined  by  equations  (3.1 O),  (3.1l)  and 

d  /•  t  ua  F(u)  du  2  ft  it  du 

^(t2  -  u2) 


h(t) 


2  d  f  t  ua  F(u)  du  2  ft  t  \ 

=  _  t" 1  —  /  _ — _  =  -  /  u  f(u)  d 

IT  dt  J  0  V(^2  -  U  *)  IT  J  Q  (t*  -  U 


(3.47) 


and  so  has  the  value  given  by  equation  (2.6).  We  therefore  obtain  the 
same  result  as  before. 

3.3.  Solution  of  Problem  (a)  in  the  case  of  constant  internal  pressure. 

We  now  consider  the  case  in  which  the  crack  is  opened  up  by  a 
constant  pressure  p  ,  The  corresponding  value  of  h(t)  is  given  by 
equation  (3.16)  so  that  if  we  make  the  substitution  (3.17)  for  g(t) 
we  find  that  the  relevant  form  of  the  integral  equation  (3 .33)  is 


4>  ( t)  -  I  <p  (u)  K(u,  t)  du  =  —  t . 
Jo  v 


(3.48) 


U2 


We  begin  by  finding  an  iterative  solution  of  this  equation  which 
is  useful  in  discussing  problems  in  which  c  >  >  1 . 

If  we  write 

00 

*(t)=-  )  hr(t)  (3.49) 

*  r^o 

then  we  have  the  iterative  scheme 

ho(t)  =  t,  hr+1(t)  =  j  hr(u)K(u,  t)du,  (r  =  0,  1,  2,  ...)  (3.50) 

i  0 

for  the  determination  of  the  functions  hr(t). 

We  now  expand  the  kernel  K(u,  t)  in  powers  of  c"1.  Using  the 
expansion 

cosh  ^  v  -  cosh  v 

_  2u  t  +  (u3t  +  ut3)  ^4  (3u5t  +  l0u3t3  +  jut5  )  v6 

c2  3  c4  180  c6 

+  lu7..^..+  7  f  +  7  U,"  t5  +  Ut7,)  v8  +  o(c-,0  ) 

2702c  8 


we  see  that  we  may  express  K(u,  t)  in  the  form 

00 

K(u,  t)  =  a  2r  +1  (u)t2r+1  , 
r=o 


where 


2u  r 

«.(»>-—[ 

t r  c  v. 


uT 


u  T 


u6T 


2  +  4  + _  6  + 


3c  60  c4  2520c 
u4  T 


] 


2u  r  u2T  u4  T  -n 

a  (u)  =  . . T  +  —  1  -  +  , , ,  { 

3  3w  c!L  4  6c2  120c4  J 


(3.51 ) 
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u  r  u  s  ~) 

a5(u)  = p6  +tt  +  •••  * 

30  v  c  4  6  c  J 

a  (u)  _  — - — t 
7  1260jt  c  4  J 


and  we  have  written 


K  (x) 

T.  =  /  -J - x  dx. 

1  1  I  (z) 


(3.52) 
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The  values  of  the  integrals  are  listed  in^Tranter  ( 1 959) . 

Prom  the  form  now  taken  by  the  function  K(u,  t)  it  follows  easily 
that  the  functions  h  (t)  may  be  written  in  the  form 


00 

hr<‘>=£  P2?  +  , 

S=1 


,2s  + 1 


(3.53) 


Suppose,  for  example,  that  c  is  sufficiently  large  to  enable  us  to  omit 
powers  of  _t  higher  than  the  seventh.  The  next  iteration  h^^ ( t)  may 
then  be  obtained  from  equation  (3.53)  where  the  new  coefficients  will  be 
determined  from  the  equation 


3  +  P^u5  +  P^uj  a.(u)du. 


P(.r+1)  a  j’jpWu  +P(3r)u: 

In  the  case  of  hQ(t),  the  initial  value,  we  take 

p  =1,  p±  =0,  1/1. 

Then,  to  find  the  coefficients  of  successive  iterations  we  need  only 
calculate  the  sixteen  quantities 


(3.54) 


(3.55) 


r  =  J  usar(u)du,  (r,  s  =  1 ,3,5,7\  (3.56) 


The  general  term  h^( t)  will  be  given  by  (3*53)  with  pj^  satisfying  the 
recurrence  relation 


p(r)_  o  P(r  -  0  ,  . 

*  "  1 ,  r  1  +  a 


p(r  -  l) 

3,  1  3 


+  a,  ,p£r"  ^  +  a  F$r  " 


(3.57) 


If  we  take 
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T  =2.50330,  T  =3.77139,  T  =23.431,  T  =302.29 

2  4  6  6 


we  find  that 

2 


r  T  T  T-n 

[2T2+-V-“t = 

4  5c  140c  7560c  J 


i  »i 


2  3 


3ir  c 


Ta  0.3382  0.0509  0.0113 


5  c  140c  7560  c  J  c  c  c 

0.0026 

+  -  y 


3»1 


2  ^2T2 

it2  c  3  L  5  21  c 


T  T  ^  0.2029  0.0364  0.0088 

— 1 - - 


540c4  27720c  6J  c3  c5  c7 


0.0022 

+  ~T"  3 


a  .-Lf 
9,1  ir2c3l 


2  r2T  T  T  T  ^  0.1449  0.0283  O.OO72 


j6J  c  3 


+  . ""■  —  +  . + 

5  _7 


7  27c2  660c4  32760c--;  c'  c'  c 

0.0019 


+  — —  , 


2  r2T  T  3T 


0.1127  0.0232  0.0061 


L  ~  T  2  3  )  +  2  +  1  *♦“  *  [  =  ,  _  +  E  "  +  4  1  + 

7,1  3irc  ^.3  11c  260 c  12600c  -J  c  c  c 

0.0016 


T  T  0.08492  0.0528  0.0243 


1  ,3 


2  5 


It  +  — s-  +  — —  1  = 

9*r  c  4  4  i0c2  280c4  J 


. .  +  —  +  . . .  , 

s  7  9  ' 

c  c  c 


2 

45  42c*  1080c4  J 


3ir  c 


0.0510  0.0377  0.0190 

— “  + - 7”  + - -  » 

c  c  o 
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a 

5,3 


2  fT4  JLu  Ta  0.0364  0.0293  0.0155 

3it2cs  t  7  54o2  1320o4J  o5  o7  o9  * 


a 

7,3 


0.0283  0.0240  0.0131 


» 


a 

1»3 


1  T  0.0264  0.0340 

-  T  +  — L-  _  - —  + - 

90ir2c7  6  10c2  c7  c9 


1  r  T  T  -v  0.0158  0.0243 


ri  x  u.ui; 

7  i  5  42c*  J  c 7 


3,5  30rr2c7  l  5  42c* 


+  9  * 

c 


1  fT  T  'J  0.0113  0.0189 


fT  T  'j  0.011 


5»s  30jt2c7  7  54c2  J  c 7  c 9 


1  rT  T  'i  0.0088  0.01547 


fls.  x  1  u.vak 

7,5  90wV  1 3  22o2J  c7  -9 


Ts  0.0081 


1,7  3780ir  o’  c' 


b6 


Te 

0.0049 

6300it2c  9 

9  * 

O 

T 

8 

0.0035 

8820 /c  9 

9  * 

C 

T 

8 

0.0027 

1 1 340  jtZc9 

9 

C 

For  instance,  it  is  easily  shown  that  up  to  terms  of  order  c 


-7 


a 

1  »5 


a  a  + 

i  *1  i  >i 


+  a  a  +  a  a 

3,1  1  ,3  3,1  1,5 


0.1144 


and  the  othera  may  be  neglected.  Since 


£ 

r,  s 


p(r)  t2s+1 
2s +1 


it  follows  from  equation  (3.19)  that 


r,  s 


2s  +  3 


and  we  can  easily  show  that 

«  (x)  =  1  +  0,3382  x3  +  0.081 5  x5  +  0.1lA4x‘+  0.0125x7  +  0(xe). 


From  the  definition  (3.22)  we  then  have 

nt(x)  =  1  -  0.3382  X3-  0.1 087 Xs-  0.1 856 x6-  0.021 9 X7  +  0(x») 
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showing  that,  in  the  notation  of  equation  (3.23),  the  critical  value  of 
the  pressure  for  a  crack  radius  a(  <  <  o)  is  given  by  the  equation 


(oo) 

^cr  ~  ^cr 


1  -  0.3382 a 3 /c3  -  0.1 087 a 3 /os  -  0.l856a6/c6  - 
-  0.0219a7/c7  +  0(a 8  /c8)j  . 


If  we  define  a  ratio  II  by  the  equation 

P  -  p 
*cr  ar 


n  = 

i 


100 


■or 


(3.58) 


(3.59) 


so  that  n  denotes  the  percentage  change  in  the  value  of  the  critical 
pressure  pcr  from  the  value  p^r  corresponding  to  a  crack  in  an  infinite 
solid,  then  if  c  >  >  a. 


ni  =  33.82a3/c3  +  10.87a!/c5  +  l8.56a6/c6  +  2.l9a7/c7  +  o(a8/c8).  (3. 59a) 

These  formulae  are  of  use  only  if  the  ratio  a/c  is  small.  When 
it  is  only  slightly  less  than  unity  the  integral  equation  (3.48)  has  to 
be  solved  numerically.  At  first  sight  the  solution  <f>  (t)  of  this 
equation  would  appear  to  depend  only  on  t  in  the  range  O  «  t  s  1  but  it 
will  be  recalled  that  the  kernel  K(u,  t),  defined  by  equations  (3.35)  and 
(3*3^) ,  depends  on  the  value  of  c,  the  radius  of  the  cylinder  expressed 
in  terms  of  the  radius  of  the  crack  as  the  unit  of  length.  The 
computations  were  carried  out  for  the  values  c  =  1.05,  1.10,  1.20,  1  .30, 
1.6667.  The  values  of  the  integral  H(x)  defined  by  equation  (3.36) 
were  first  calculated  vising  Weddle's  rule  for  numerical  integration  and 
the  integral  equation  itself  was  solved  using  the  method  of  Pox  and 
G-oodwin  (l  95  ).  The  calculations  were  carried  out  on  the  DEUCE  computer 
in  the  Computing  Laboratory  of  the  University  of  Glasgow. 

The  results  of  the  calculations  of  $(t)  are  shown  in  Table  3  and 
graphically  in  Pig.  7.  It  will  be  observed  that  even  for  moderately 
small  values  of  c  (for  instance  greater  than  1 .20)  the  graph  of  the 
function  <#>(t)  differs  very  little  from  a  straight  line. 

Using  these  values  of  <p  (t)  and  the  formulae  (3.18)  and  (3.19)  we 


c-1.05 


c-1.10 

c-1.20 
c«1.30 
c- 1.667 
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ratio  W/W  where  W 

oo  oo 

,  i.e.  to  a  value  0  of 


Table  4  and  graphically  in 

Fig*  8* 

To illustrate  the  variation  of  the  critical  pressure,  p„r,  necessary 
to  cause  the  spread  of  the  crack  when  its  radius  is  a  the  ratio  n 
defined  by  equations  (3«59),  (3.2 3),  (3.22)  and  (3.19)  was  calculated. 

The  results  are  shown  in  Table  4  and  graphically  in  Fig,  9.  From  this 
diagram  it  follows  that  the  drop  in  the  value  of  the  critical  pressure 

(oo\ 

Pcr  from  the  value  pQr  corresponding  to  a  crack  in  an  infinite  solid  is 
less  than  5$  if  c  >  2a  and  is  less  than  10$  if  c  >  1 ,4a. 

If  we  look  at  the  variation  with  a/c  of  the  stress  intensity  factor 
N  we  get  much  the  same  kind  of  picture.  To  illustrate  the  variation  of  N 

we  calculated  n  the  percentage  increase  in  the  stress  intensity  factor 
due  to  the  effect  of  the  finite  radius  of  the  cylinder.  The  results  are 
shown  in  the  last  row  of  Table  4  and  graphically  in  Fig.  10.  The 
effect  of  the  finite  value  of  the  radius  is  more  pronounced  in  this  case. 
For  instance  if  c  <  1.2a,  the  change  in  the  stress  intensity  factor  is 
greater  than  40$.  However  it  is  less  than  5%  if  c  >  2a  and  is  less  than 
10$  if  c  >  1 .7a. 

It  is  of  interest  to  compare  these  results  with  those  obtained  in 
the  analogous  problem  in  plane  strain.  If  we  plot  the  ratio  n  i  as  a 
function  of  x,  the  ratio  of  the  area  of  the  crack  to  the  cross-sectional 
area  of  the  cylinder  (i.e.  x  =  a2/ c 2)  we  get  the  curve  shown  in  Fig.  11. 

The  curve  n  shown  in  the  same  diagram  shows  the  value  of  the  percentage 
change  in  p  in  the  plane  strain  problem  but  here  x  is  taken  to  be  the 
ratio  of  the  length  of  the  crack  to  the  width  of  the  strip.  In  Fig.  12 

we  show  a  comparison  of  the  percentage  increase  in  the  stress  intensity- 
factor  in  the  axisymmetric  case  with  n  the  corresponding  quantity  in  the 
plane  strain  case,  x  being  given  the  same  interpretation  as  before.  Using 
this  basis  of  comparison  we  see  that  as  far  as  an  "engineering"  approx¬ 
imation  is  concerned  the  size  effect  in  the  axisymmetric  case  can  be 
gauged  from  that  in  the  plane  strain  case. 


easily  calculate  the  variation  with  a/c  of  the 
denotes  the  value  of  W  corresponding  to  c  =  « 
the  ratio  a/c,  so  that 

8  P*0  -  /)a3 

W  =  - - -  . 

“  3  E 

The  results  of  these  calculations  are  shown  in 


Figure  10 


Fig.  10  The  variation  «ith  c/a  of  ^  the  percentage  increase  in  the 
stress  intensity  factor  due  to  the  effect  of  finite  radius  of  the 
cylinder. 


Figure  11 


Fig.  11  A  comparison  of  the  axi symmetric  solution  and  the  plane  strain 
solution,  rr^  is  the  percentage  change  in  the  value  of  the  critical 
pressure  in  the  axisymmetric  case  with  jc  equal  to  the  ratio  of  the  area 
of  the  crack  to  the  cross-section  area  of  the  cylinder 5  tt  is  the  analogous 
quantity  in  the  plane  strain  problem  with  jc  equal  to  the  ratio  of  the 
length  of  the  crack  to  the  width  of  the  strip. 
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Figure  12 

Fig.  12  A  comparison  of  the  axisymmetric  solution  and  the  plane  strain 
solution,  n  is  the  percentage  increase  in  the  stress  intensity  factor 
inrthe  axisymmetric  case,  n  that  in  the  plane  strain  case,  (x  is  defined 
as  in  Fig.  ll). 
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3.4..  Solution  of  Problem  (b). 

In  the  second  of  the  two  problems  we  shall  consider  the  conditions 
on  the  crack  surface  are  the  same  as  in  the  first  problem  but  the 
conditions  on  the  curved  surface  of  the  cylinder  are  different.  Instead 
of  assuming  that  the  normal  component  of  the  displacement  vanishes  we 
assume  that  the  normal  component  of  the  stress  vanishes,  i.e.  the 
boundary  conditions  (3.4)  are  replaced  by  the  conditions 

a  -  o  =0.  p  =  c.O$z<oo  (3»6o) 

pp  p  z  *  *  ’ 


If  we  take  the  form  (3.6)  for  the  stress-function  X  (  P  ,  z)  we 
find  that  on  the  surface  p  =  a  the  stress  component  0  assumes  the  form 

jjO  =  ^  I  ji(5)jo«lo(o5)  -  1,(05)]  +b(5)[(3  -  2n)cslo(c£)  - 
-  4(1  -  T)  )l,(c€)  -  c2  €2  1,(0^)]  J  oos(€  z)d€ 

,  r  J  (fo) 

-  2p  /  F(e)  e“  5z  J  (  C  c)  -  (1  -  2  TJ  )  -J - Jz  x 

Jo  C  0  {c 

[JoUo)_L^]]aj  . 


J.(s  0) 

S' 

If  we  take  the  Fourier  cosine  transform  of  both  sides  of  this  equation  and 
make  use  of  the  fact  that  E  pp]p  ac  =  0  and  hence  that  ^olD^c 


]' 


£  we  find  that 


A(5)  £o5I0(c5)  -l,(c5)J  +b(5)[(3  -2n)o£I0(c5)  -4(1  -n)  x 
I,(o€)  -  c2  =  c5i3  +  rji,-ig  (3.61) 

where  i,  and  ig  are  defined  by  equations  (3.28)  and  ig  is  defined  by  the 
equation 

45  sp(  s)  j0(  S°)  dS 

. -  T“  • 


= 


(5  ♦  S  ) 


(3*62) 
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Substituting  the  expression  (3.9)  for  P(  )  into  equation  (3.62), 
interchanging  the  order  of  the  integrations  and  making  use  of  well-known 
integrals  involving  Bessel  functions  (Cf.  pp.10,  11  of  vol.ii  of 
Erdelyi  (l  954)  ),  we  find  that 

i3  =  —  j  g(u)  j  sinh(  gu)  l~K  o(c  5)  -cgK^c?)"!  +  5ucosh(£u)x 
rr  <*  o  L  *“  1  -J 

Ko(c?)J  du.  (3.63) 

Solving  the  equations  (3.6l),  (3.27)  we  find  that  the  boundary 
conditions  ( 3  •  ^  0)  are  satisfied  if  we  choose  A(  %) ,  B(  g  )  to  be  given  by 
the  expressions 

?irA(5)&(cg)  =[h(o€)  -  2  +  2nJ  T  g(u)  gu  cosh(  gu)  du 


+  ["(3  -  2  v  )H(c  5)  -  4+  47?  -  cZg^j  j  g(u)sinh(  £u)  du, 

(3.64) 

£  jtB(  g  )Cr(c  £)  =  J  g(u)  g  u  cosh(  £  u)  d  u  +  jj  -  H(c  $  )j  j"  g(u)  sinh  (  5  u)  du 

(3.65) 


where  the  function  G  is  defined  by  the  equation 

G(s)  =  s2I2(s)  -  (2  -  2 77  +  sz)l2(s)  (3.66) 

‘0  1 

and  the  function  H  by  the  equation 

* 

H(s)  =  ■*!  (s)K  (s)  +  (2  -  277  +  s2)l  (s)K  (s).  (3.67) 

00  11 

If  we  substitute  from  equations  (3.64)  and  (3.65)  into  equation  (3.11 ) 
we  find  that  this  last  equation  reduces  to  the  integral  equation 

g(t)  +—  f  L(t,  u)g(u)du  =  h(t),  Ok  t  <1,  (3.68) 

v  Jo 


for  the  determination  of  the  function  g(t),  the  kernel  L(  t,  u)  being 
defined  by  the  equation 


58 


L(t,  u)  g(c  g^O^0  s)  -  i]  oosh(g  u)slJih(  gt)  +  gtcosh(gt)x 

sinh(i-u)[]  +  [2H(cg)  -3+2  rj  -  o2  g^]sinh(  £u)sinh(  gt) 

-  £Z  utcosh(gu)  cosh(gt)Jd£.  (3.69) 

The  equations  (3.68)  and  (3.69)  are  identical  with  those  derived  by 
Oollins  Q)f.  equations  (3. 13)  and  (3.14)  of  Collins  (1962)^}. 

For  small  values  of  the  constant  c"1  it  is  possible  to  solve 
the  integral  equation  (3.68)  by  an  iterative  procedure.  To  derive  this 
solution  we  need  first  to  derive  an  expansion  of  L(t,  u)  in  powers  of  c " 1 . 
We  find  that 


L(t 


00  A 

,  u)=  V  c-2r"3  y  D  t2n+1  2r-2n+1 

*  '  L-i  L-i  r-  n,n 


(3.70) 


r=o 


n=o 


where 


B„,  n  - =  20"  +  n  *  2)B.,n-  [2<“  *  ’>(»  *  h  .  1  -1]^  - 

-  (.  ♦  l)(2.  ♦  (3.71) 


with 


S  =  E  = - - - [ 

m,n  n,m  n  (2m  +  1 )l(2n  +l)N( 


»  2m+2n+2  T 


s  H(s)d  1 

0  G-(s) 


(3.72) 


and 


F  =  F 
m,n  n,m 


w(2m 


- 2 - f 

+  l)i(2n  +  l)I J i 


«  2m+2n+2  , 

s _ as 


<X») 


(3,73) 


For  a  prescribed  value  of  the  Poisson  ratio  rj  the  coefficients  E _ and 

m,n 

Fm  n  have  to  be  evaluated  numerically  for  various  values  of  m  and  n  depending 
on  the  order  of  the  terms  retained  in  the  expansion  on  the  right-hand  side 
of  the  equation  (3.70). 
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3.5.  Solution  of  Problem  (b)  in  the  Oaae  of  Constant  Internal  Pressure. 

We  now  oonsider  in  more  detail  the  oaae  in  which  the  internal 
pressure  in  the  crack  is  p  so  that  f(p)  =  pQ  /2p  and  h(t)  =  Po 
and  the  integral  equation  (3.68)  becomes 


2  fi  p  t 

g(t)  +  -  /  L(t,  u)g(u)du  =  — ,  0  «  t  «  1 . 

v  **  o  it  ia 

Collins  has  obtained  the  approximate  solution 


(3.74) 


(3.75) 


by  the  method  of  iteration. 

Now  by  equations  (3.1 3)  >  (3.14)  we  find  that  the  energy  of  the  crack 
is  given  by  the  equation 

W2  =  4jrpo(l  -  77  )  J  tg(t)dt  (3.76) 

*  0 

so  that  in  our  units  (in  which  the  radius  of  the  crack  is  taken  to  be  the 
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unit  of  length) 
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In  conventional  units  this  becomes 


W  = 
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+  oCaVc’1)].  (3.77) 


Applying  the  Griffith  criterion  (3.20)  we  find  that  the  critical  value  of 
po  is  given  by  the  formula 


In  the  case  77  =  0,3  we  find  that  this  formula  takes  the  farm 


pcr  =  PS  £  "  °‘2626  a’/c3  -  0.0340  a5/ c5  -  0.4570  a7/c7  +  0.0176ae/c8 

,0(aV«’)]  (3.79) 

If  we  define  a  percentage  change  II  in  the  value  of  the  critical  pressure 

§  1- ■  V  ^ 

Por  from  the  value  p^”in  the  infinite  case  by  equation  (3*59)  then  we 
find  that 

n2  =  26.2 6  a3/c3  +3.40  a*/o5  +  45.70  a7/c7  +  0(a8/c8).  (3.80) 

Formulae  (3.79)  and  (3.80)  should  be  compared  with  their  counterparts 
(3.58)  and  (3.59)  derived  for  the  boundary  conditions  of  Problem  (a). 

When  a/c  is  not  very  small  we  again  have  recourse  to  numerical 
integration  of  the  integral  equation.  To  facilitate  the  numerical 
calculations  we  transform  the  integral  equation  (3.68)  by  making  the 


substitutions 


6l 


Po° 


2  o 


g(t)=—  K(t,  u)  = —  L(ot,  cu),  a  =  c 

ff/u  rr 

when  it  reduces  in  the  case  h(t)  =  p o^/(  n  n)  to  the  form 
4>  (t)  +  f  K(t,  u)  <f>  (u)  du  =  t,  0  s  a  <  1 

J  0 


(3.81 ) 


(3.82) 


where  the  kernel  is  now  given  by  the  equation 

K(u,  t)  =  jj2H(x)  -3+2  r?  -  x^  sinh(u  x)sinh(t  x)  - 

v 

-  ut/  cosh(u  x)cosh(  t  x)  +  xQh(x)  -  ac»  cosh(u  x)sinh(t  x)  + 

+  t  sinh(u  x)cosh(t  x  0]sfe  (3.83) 

with  H(x)  and  G-(x)  defined  by  equations  (3.66)  and  (3.67). 

Prom  equations  (3.13)  and  (3.1 4-)  we  find  that  the  energy  W 

2 

required  to  open  the  crack  is  given  (in  conventional  units)  by  the 
equation 


W  = 
2 


8  p2n(l  -  O  _  fa/  c 


E 


/  t  0  (t)  d  t 

i  O 


(3.84) 


where  a  is  the  radius  of  the  crack,  so  that  the  Griffith  criterion  (3.20) 
that  the  crack  may  spread  leads  to  the  formula 


with 


p  =  p(^  Q  ,  Wo)  . 

*cr  ■'or  2  '  * 

(3.85) 

Hg(a)  =  a/  <f>  (a)J  ^ 

(3.86) 

for  the  critical  value  of  the  applied  pressure,  p^  being  defined  as  before. 

The  percentage  change  IIg  in  the  value  of  p^J  is  therefbre  given 

by  the  formula  r 

nt=  ioo  ji  -  n2 (a/o)j 


(3.87) 
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In  a  similar  way  equation  (3.12)  leads  to  the  expression 
4(l  -  n2  )p  a 

(p)  = - —  ?a  <*),  0.  x  <  a 

jrB 


with 


a 


e/c,  x  = 


P/°»  F  a  (x) 


1  rot  $(t)dt 

a  L  A/( t2  -  x2) 


(3.88) 


(3.89) 


for  the  normal  displacement  of  the  surface  of  the  crack. 

4(l  -  n  *)po  a 
b  B  1  ■  111 

irE 


If  we  write 

(3.90) 


which  is  the  value  of  u  at  the  origin  of  coordinates  in  the  case 
c  »  a  we  see  that  we  can  write  equation  (3.88)  in  the  form 


u°  (p)  =  b  P  (x).  (3.91) 

i 

The  percentage  increase  n  g  in  the  stress  intensity  factor  can  be 
calculated  by  means  of  equation  (3.60) .  We  find  that 

n2  -  1  J  x  100,  a  =  a/c.  (3.92) 

The  various  functions  were  computed  for  the  value  77  =0.25. 

The  variation  of  <p  (t)  with  t  for  six  values  of  a  is  shown 

in  Table  5  and  graphically  in  Pig.  13*  while  that  of  P  (x)  for  0  $  x  <  1 

and  the  same  values  of  a  is  shown  in  Table  6  and  Pig.  14.  The 

variation  of  the  function  Q  2  (a/ c)  i.e.  of  the  ratio  VOI/v£/!  is  given 

by  the  second  row  of  Table  7.  Pig.  15  shows  the  variation  of  the 

critical  pressure  pQr  with  the  ratio  a/a.  The  values  of  the  other 

quantities  of  physical  interest  Ha  and  n ^  are  also  listed  in  Table  7. 

To  illustrate  the  differing  effects  of  the  two  kinds  of  boundary 

conditions  on  the  cylindrical  surface  -  those  considered  here  and  those 

considered  in  problem  (a)  -  the  variation  of  II  and  n  is  shown  in 

“  12 

Pig.  16  and  that  of  n  and  n  is  shown  in  Pig.  17. 

1  2 

Using  either  of  these  quantities  as  a  measure  of  the  effect  of  the 
finite  radius  of  the  cylinder  we  see  that  if  a  >  |  o  the  effect  is  more 
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Figure  ll* 

Fig.  1U  The  variation  withx  and  a  of  the  function  Fjx)  for  constant 

p  and  tq  =  0.25. 


0  0.2  0.4  0.6  OB  t.0 

Figure  16 

Fig.  16  The  variation  with  a/c  of  the  percentage  increase  in  the  critical 
pressure  for  constant  pressure  on  the  crack  surfaces.  The  curve  vr^  corresponds 
to  the  case  in  which  up  =  0  on  the  cylindrical  surface  and  is  independent  of 
t)j  the  curve  tt^  corresponds  to  the  case  in  which  app  =  0  on  the  cylindrical 
surface  andr)=  0.25. 


Figure  17 

Fig.  1?  The  variation  with  a/c  of  the  percentage  incre  ase  in  the 
stress  intensity  factor  for  constant  pressure  on  the  cra_ck  surfaces. 

The  curve  n^  corresponds  to  the  case  in  which  up  -  0  on  the  cylindrical  sur 
face  and  is  independent  of  the  curve  ng  corresponds  t—o  the  case  in 
which  <7  =  0  on  the  cylindrical  surface  and  t)  ■  0.25. 


pronounced  in  the  case  in  which  the  radial  component  of  the  surface 
displacement  of  the  cylinder  is  zero. 
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3.6.  Solution  of  Problem  (b)  in  the  Oase  of  a  Variable  Internal  Pressure. 

To  illustrate  the  effect  of  the  application  of  a  variable  pressure  to 
the  surfaces  of  the  crack  the  numerical  calculations  were  repeated  for  the 
case  in  which  2  . 

f(p)  - *  )  ,  0  s  p  <  1 .  (3.93) 

2p  \  c  J 


In  this  case  the  free  term  of  the  integral  equation  (3.68)  assumes  the 
form 


v\i 


(3.94) 


If  we  insert  this  value  for  h(t)  into  equation  (3.68)  and  make  the 
transformations  (3.81 )  we  find  that  the  function  <+>  ( t)  is  a  solution  of 
the  integral  equation 


K(t,  u)  <f>  (u)du  =  t  -  j  t3 


(3.95) 


where  the  kernel  K(  t,  u)  is  given  by  the  equation  (3.83). 

Numerical  solutions  of  this  equation  were  obtained  for  the  four 
values  0.2,  0.5,  0.7  and  0.9  of  a  =  a/o  and  for  ti  =  0.25.  The  results 
are  given  in  Table  8  and  are  illustrated  graphically  in  Pig.  18.  The 
function  PQ  (x)  defined  by  equation  (3.89)  was  also  calculated  for  these 
four  values  of  a  .  The  results  are  shown  in  Table  9  and  graphically  in 
Pig.  19.  Here  the  upper  curves  show  the  distribution  of  pressure  on  the 
crack  surfaoe  and  the  lower  curves  give  the  resulting  shape  of  the  crack; 
the  length  b,  the  depth  of  the  crack  in  the  case  a  «  c,  is  defined  by 
equation  (3 • 90) . 

The  quantities  0g  and  ng  which  are  defined  by  equations  (3.86) 
and.  (3.92)  and  from  which  can  be  derived  the  aritical  pressure  and  the 
percentage  increase  in  the  stress  intensity  factor,  were  calculated  far 
these  four  values  of  the  ratio  a/o.  The  results  are  shown  in  Table  10. 

In  comparing  these  results  with  those  corresponding  to  oons tan t  internal 
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Figure  18 


Fig.  18 
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The  variation  with  t-  and 

p  2 

p  (1-  p  /c  )  on  the  crack 
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surfaces  and  r)  =  0.25* 


Figure  19 

Fig.  19  The  shape  of  the  crack  for  four  values  of  a/c  when  the  pressure 
is  p  (l-p2/c2)  and  T)  =  0.25.  The  upper  curves  show  the  distribution  of 
pressure  on  the  crack  surface  and  the  lower  curves  the  corresponding  shape 
of  the  crack.  [b  -  U(l~n2)poa/ffE). 
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pressure  (Table  7  above)  it  should  be  noted  that  in  the  second  case  the 
average  pressure  on  the  crack  is 

2 p 0  a" 2  I  p  (l  -  p  2/o  2  )dp  =  pQ  fl  -  £(a/c)  H 

J  o  J 

so  that  if  a  is  nearly  equal  to  c  the  average  pressure  is  appreciably 
less  than  p  . 

*  n 
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